Optical signatures of the tunable band gap and valley-spin coupling in silicene 
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We investigate the optical response of the silicene and similar materials, such as germanene, in 
the presence of an electrically tunable band gap for variable doping. The interplay of spin orbit 
coupling, due to the buckled structure of these materials, and a perpendicular electric field gives 
rise to a rich variety of phases: a topological or quantum spin Hall insulator, a valley-spin-polarized 
metal and a band insulator. We show that the dynamical conductivity should reveal signatures 
of these different phases which would allow for their identification along with the determination of 
parameters such as the spin orbit energy gap. We find an interesting feature where the electric 
field tuning of the band gap might be used to switch on and off the Drude intraband response. 
Furthermore, in the presence of spin-valley coupling, the response to circularly polarized light as 
a function of frequency and electric field tuning of the band gap is examined. Using right- and 
left-handed circular polarization it is possible to select a particular combination of spin and valley 
index. The frequency for this effect can be varied by tuning the band gap. 
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I. INTRODUCTION 



The experimental realization of two-dimensional crys- 
tals has sparked a flurry of research activity. The proto- 
typical example was the isolation of graphene j 1 ^ a single 
layer of carbon atoms, with its remarkable physics due to 
the Dirac fermionic nature of its charge carriers. Much 
research has ensued^— and part of the focus is on the 
potential to manipulate the physics of these charge carri- 
ers to benefit technological applications. In graphene the 
charge carriers show a linear dispersion near the Fermi 
leve£ with chiral properties attributed to the topological 
aspect of two sublattices in a hexagonal crystal lattice. 
This leads to a mapping of the low energy nearest neigh- 
bor tight-binding Hamiltonian on to a Dirac equation for 
massless fermions which exist at K and K' = —K in 
the hexagonal Brillouin zone£ The conical disperison at 
these K points are referred to as valleys and many prop- 
erties in graphene are spin degenerate and valley degen- 
erate. As the spin-orbit interaction in graphene is quite 
small^ and therefore neglected, interesting effects, such as 
inducing spin-valley coupling discussed here, will not be 
evident. This is unfortunate as much interesting physics 
and potential for spin manipulation will not be manifest 
in graphene. Given some of the limitations of graphene 
in this regard, attention has turned to other 2D crys- 
tals which have graphene-like features but stronger spin- 
orbit coupling. One such material is silicene, a monolayer 
of silicon atoms arranged in a honeycomb lattice, as in 
graphene. As silicon is already used in electronics and a 
large silicon manufacturing industry is in place, it bears 
closer inspection. Silicene has only recently been made in 
ribbons 1 ^— and an experimental literature has yet to be 
developed. However, numerical calculations have shown 
that it should be possible to make silicene suitable for 
transfer to a substrate where it might be gatedi 13 i 14 As 
silicon has a stronger spin-orbit coupling than graphene, 



it will exhibit more clearly an energy gap in the band 
structure. It has been argued that in this circumstance, 
the Dirac nature of the fermions makes a topological insu- 
lator (TI) . One of the interesting predictions for silicene 
involves that of inducing a tunable band which 
would give rise to a transition between a topological in- 
sulator and a band insulator (BI) . 14 ' 15 This comes about 
due to the fact that in silicene, the size and distance 
between silicon atoms leads to a buckling of the struc- 
ture where one sublattice is shifted upwards out of the 
2D plane relative to the other (see Fig. [TJ. This means 
that an applied electric field perpendicular to this plane 
changes the sublattice potential and produces a tuning 
of the band gap in the dispersion^ Furthermore, the 
bands are spin-split in a valley-dependent way allowing 
for an examination of spin- valley coupling. 16 ^ Other exam- 
ples of systems discussed in the literature which are pre- 
dicted to show spin-valley coupling are M0S2 and other 
group- VI dichalcogenide o 17 i 18 and a monolayer of germa- 
nium (germanene) 1 ^. Currently, M0S2 has been the sub- 
ject of several experimental investigations, for example 
Ref i 18 ' 20 " — , while silicene is still in its infancy and ger- 
manene has yet to be synthesized. Nonetheless, we will 
use the examples of silicene and germanene as a models to 
examine issues associated with spin-valley coupling and 
optical signatures of a TI to BI transition. 

In this work, we study the finite frequency conductiv- 
ity of a material such as silicene which can be induced to 
have spin-valley coupling and a tunable band gap. We 
find that there are characteristics in the far infrared or 
THz conductivity which would allow for the differentia- 
tion of the TI versus the BI. With tuning of the band gap, 
these features could be used to identify the critical point 
of transition and allow for an experimental determination 
of the strength of spin orbit coupling in these systems. 
We also consider the effect of finite doping and circularly 
polarized light as possible probes of spin- valley coupling. 
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FIG. 1. (Color online) The crystal structure of silicene is 
based on the honeycomb lattice shown here, however, due 
to the size of the silicon atom, the the honeycomb lattice is 
buckled with the A and B sublattices shifted out of a com- 
mon two-dimensional plane. [Crystal structure plotted using 
VESTA^.] 



Previous issues associated with circular dichroism are re- 
visited in view of performing broadband spectroscopy. 

Our paper is organized in the following manner. In 
Section II, we briefly describe the basic Hamiltonian and 
approximations used in this work and outline the nature 
of our calculation of finite frequency conductivity. Fol- 
lowing this, in Section III, we feature our results for the 
absorptive part of the longitudinal optical conductivity 
and provide a discussion of the physics. In Section IV, 
we discuss the case of circularly polarized light which il- 
lustrates the potential for spin- valley polarization and its 
the tuning by electric field. Our conclusions are summa- 
rized in Section V. 



II. THEORETICAL BACKGROUND 

Recently, there have been theoretical works which have 
calculated the band structure of silicen o 14 ' 19 using den- 
sity functional theory including effects due to spin-orbit 
coupling (SOC) and a perpendicular electric field. The 
spin-orbit band gap has been found to be about 1.5 meV 
and this can be increased to 2.9 meV under strain^ This 
latter study has also predicted a SOC band gap of 23.9 
meV for 2D germanium or germanene. In the work of 
Drummond et ali 4 ., a SOC gap of 1.4-1.5 meV is also 
found for silicene and with a perpendicular electric field, 
the resulting gap may be tuned up to tens of meV before 
other transformations occur. This work along with that 
of Ezaw a 15 i 16 ' 25 suggest that much of this behavior may 
be captured by a simple low energy Hamiltonian written 




FIG. 2. (Color online) The presence of spin orbit coupling 
and a perpendicular electric field, parametrized by Aso and 
A 2 = E z d, respectively, gives rise to spin-split bands about 
the K point, with two gaps as shown in the upper left, one 
of which may be tuned to zero at A z = Aso as shown in the 
lower left. The bands at K' are reversed from those at K 
(right frames). 

about the K point labeled 

H K( = v(£,k x r x + k y T y ) - t-AsoZzh + \& z t z , (1) 

where £ = ± is the valley label for the two K-points, 
K± which were previously referred to as K and K' . The 
first term is the Hamiltonian for Dirac electrons with 
velocity v (we take K = 1), the second term is the Kanc- 
Mele SOC term with SOC gap Ago and the last term is 
associated with the perpendicular electric field E z which, 
when applied across the buckled structure, gives rise to 
an A-B sublattice asymmetry with on-site potential A 2 = 
E z d where d is the perpendicular distance between the 
two sublattice planes. The matrices fj and &i are the 
Pauli matrices. The former act on the pseudospin space 
associated with the A and B sublattices and the latter 
are for real spin with cr z , the z-component. In the papers 
by Ezawa, a Rashba SOC is included in addition to the 
intrinsic SOC used here. As the Rashba SOC is an order 
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of 10 smaller in magnitude, we neglect it. Therefore, the 
full 8x8 matrix spanning both K points is block diagonal 
in 2x2 matrices labeled by valley index £ = ± for K and 
K' and spin index a = ± for spin up and spin down, 
respectively: 



H, 



v(£k x + ik y ) 



so 



iky) 



(2) 



The eigenvalues are evaluated from the Hamiltonian 
above to be: 



± 



v 2 k 2 + i(A a 



^Aso) 2 - 



(3) 



These are plotted in Fig. [2] for a value of A 2 = 0.5Aso 
where the bands are spin-split, reversed at the other K 
point, and represent a topological insulator. Also shown 
is the case where A 2 = Ago- At this critical point, 
the gap of one of the spin-split bands closes to give a 
Dirac point while at the other K point it is gapped and 
it is the other spin band which has no gap. This has 
been termed a valley-spin-polarized metal (VSPM)/2£ For 
A 2 > Aso, the spectrum becomes fully gapped again but 
the system is a band insulator albeit with unusual chiral 
properties^ In the case of A 2 = 0, the bands are spin 
degenerate with a single band gap of Aso- In Fig. [3j we 
show the total density of electronic states for the bands 
shown in Fig. [5J The density of states reflects the two 
gaps that open in the system depending on the strength 
of E z . At the critical value where A 2 = Asoj which is 
the VSPM state, the low energy density of states is lin- 
ear out of zero energy and jumps at a value equal to the 
sum of A 2 + Aso to further linear behavior with double 
the slope (upper righthand frame of Fig. [3]). This jump 
reflects the spin-split bands. For finite A 2 ^ Aso, a low 
energy gap appears at a value of the difference of these 
two energies and a jump at higher energy remains (upper 
lefthand frame). If E z = 0, there is only the spin-orbit 
gap and the spin degeneracy in this case gives just the one 
gap edge with no further jumps in the density of states. 
This behavior is embodied in the following formula for 
the density of states: 



N{e) = Y J N ai {e 



where 



2ttH 2 v 2 



e(2|e| - |A ct€ | 



(4) 



(5) 



and A(j£ = A 2 — er£ Aso j with a = ± indexing the spin 
up/down bands at the Kg point. A plot of the energy 
gaps for the spin-resolved bands are shown in the lower 
frame of Fig. [3J where the solid curve follows the notation 
of Fig. [2] and indicates the gap for the spin down band 
and the dashed line is for the spin up band at K. The 
solid and dashed lines would reverse in the plot for the K' 
point, as seen in Fig. [5J The region of electric field that 




FIG. 3. (Color online) Upper frames: The total density of 
states for the band structures illustrated in Fig. [2] The right 
frame is the case for the VSPM state. The variation with E z 
of the two edges shown in the left frame will follow the trend 
shown in the lower frame. Lower frame: The evolution of gap 
in the spin down (solid curve) and spin up (dashed curve) 
band about the K point as a function of A z = E z d. The 
regions of band insulator (BI) versus the topological insulator 
(TI) are indicated by shading and the dotted line marks the 
valley-spin-polarized metal (VSPM) state. 



gives |A 2 | < Aso has been shown to pertain to a TI of 
the Kane-Mele type (or a quantum Hall spin insulator) 
and for |A 2 | > Aso, the system is a BI. The critical 
point at which \A Z \ = Aso is the VSPM state discussed 
previously. It is clear from this figure that in the BI 
state, both gap edges in the DOS will increase with the 
magnitude of the perpendicular electric field, whereas in 
the TI case, the higher one will increase with electric field 
while the lower one decreases, which would be a method 
of differentiating between the density of states of the BI 
versus the TI. 

Such considerations also will apply to the dynami- 
cal conductivity which we evaluate by the usual Kubo 
formalism.— ~— The zero temperature, absorptive part 
of the conductivity at photon frequency Q is given from 
an evaluation of 



e 

2ft 



did 



fi\ — n 



d 2 k 



2tt J (2tt) : 



Tr 



where fj, is the chemical potential, 



ViA(k, uj + £l)ijA(k, ui) 
(6) 

i>i = dH/dki, with 
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i,j referring to spatial coordinates in 2D, and A(k,uj) is 
the spectral function of the Green's function found from 
G~ 1 (z) — zl — H through the relation 



duj Aij (oj) 
2tt z — uj 



(7) 



Using Eq. ([B]) and the appropriate Vi operator, we obtain 
the expression for the real part of the zero temperature 
longitudinal conductivity 



Rea xx (0) = ^Re<!(fi), 



(8) 



where the response associated with a particular spin in a 
particular valley is given as 



Re<|(0) 



e 

20 



(2tt) 



d 2 k 



(9)" 



with 



Ali(k,uj) = 7rJ2 1 

a=± ^ 

^(*, W ) = 7rWl 



vA, 



2|W 



a=± 



2k 



<5(w + a|e^|), (10) 



S(io + a\e^\), (11) 



and £3-5 is given by Eq. ([3]). From this we are able to 
obtain analytically, the real part of the longitudinal con- 
ductivity and its Kramers-Kronig-related imaginary part. 
These are given as 

ReoZKn) =a | 4/X 4 ~^^ ffi)6(2| M l - |A^|) 



Q 

2 



27rOO r 



1 

47T 



6(0 -fi c )}-, (12) 
/(«) 



1 



(13) 



where fl c = max(2|/i|, |A^|), /(fi) = ln[|0 + Q c [/|f2 - 
O c |], and (To — e 2 /Ah. Likewise, the imaginary part of 
the zero temperature transverse or dynamical Hall con- 
ductivity for a given spin and valley is given by 

Im<r4(0) = — / — 



yy ' 20./i,.i_o 2vr J (2tt) 2 



A?f (fe, w)^||(fc,w + SI) - A*(k,u)A'*(k i u> + n) 
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FIG. 4. (Color online) The real part of the dynamical conduc- 
tivity Kea xx (Sl)/ao for fi = and varying (a) A z /Aso < 1 
and (b) A z /Aso > 1. Inset: A plot of the minimum and 
maximum band gaps, reproduced from Fig. showing the 
evolution of the absorption peaks with A z . 



where A\\{k,u)) and A 2 2(k,uj) are given by Eqs. (jTUJ) 
and (TTTT) , respectively. This can be solved analytically 
for the imaginary part and the Kramers-Kronig-related 
real part. These are 



Im<|(0) = -e<7 Q -^e(0-0 c ) 



A 



Re<f(n)=fr -£/(n) 



20 



27rO" 



(15) 
(16) 



Similar analytic forms have been found in other works 
for other systems i 29 ' 30 With these formulas in hand, we 
can now discuss the longitudinal conductivity and the 
conductivity for circular polarization associated with the 
different valley and spin degrees of freedom. Our aim is 
to identify signatures of of the VSPM state and the TI 
to BI transition. 



III. LONGITUDINAL OPTICAL RESPONSE 

In this section, we discuss the results for the real 
part of the frequency-dependent longitudinal conductiv- 
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FIG. 5. (Color online) The real part of the dynamical con- 
ductivity Kecr xx (fl) /cto for /u/Aso = 0.25 and varying (a) 
A 2 /A S o < 1 and (b) A z /Aso > 1. Inset: A plot of the 
minimum and maximum band gaps, reproduced from Fig. [3j 
showing the evolution of the absorption peaks with A z but 
with a lower bound of 2/i. 

ity. For our calculations shown in the figures, we have 
chosen to evaluate Eqs. © and ©, where we have re- 
placed the delta functions in the spectral functions with 
Lorentzians according to 6(x) — > (ti/it)/(x 2 + rj 2 ) with 
rj, a broadening parameter taken to be rj = O.OlAgo- 
(Note that the Kramers-Kronig related piece of the delta 
function in Imcr^l (f2) will have to be replaced with that 
for the Lorentzian form.)^ This mimics impurity scat- 
tering in the system, with the transport scattering rate 
l/^imp = 277. This is done to provide more realistic 
looking curves, however, the analytical formulas above 
work very well if the delta function is replaced with a 
Lorentzian with halfwidth associated with the transport 
impurity scattering rate, instead of 77 which is the quasi- 
particle scattering rate in the spectral function ; 26 ' 28 

For the case of charge neutrality, where fi = 0, we have 
results shown in Fig 2] The conductivity is entirely based 
on interband transitions. For A z = 0, we have only spin 
orbit coupling, the bands are degenerate in spin with a 
single gap, and there is one absorption edge at Ago which 
shows a peak before decaying to the universal background 



expected for a monolayer graphene-like system. With fi- 
nite A z < Ago (the TI regime) , two edges appear due to 
the spin-split bands giving rise to two gaps as discussed 
in Fig. [3J (Note that transitions can only occur between 
bands with the same spin index.) The absorption edges 
track the gaps shown in Fig. [3] bottom frame (and re- 
produced here in the inset for positive A z ) such that the 
edges move oppositely in frequency - a signature feature 
of the TI region. At the VSPM point (A z /A so = 1), the 
VSPM gives the famous Dirac massless fermion conduc- 
tivity, which is flat, but the amplitude is only one half 
that of cto which is different from the case of graphene 
where it is cto . At the energy scale of 2 Ago there is now 
an absorption edge leading to a peak structure, as before, 
with the onset of the interband transitions between the 
gapped bands. This curve is a signature for the VSPM 
state. For the region for the BI (A z > Ago); two inter- 
band edges occur, the peaks of which both move outward 
according to the sum and difference of A z and Ago, but 
with a constant separation in energy of 2 Ago- The be- 
havior of the absorption edge peaks from the TI to BI 
regimes tracks the arrows of the inset, i.e., as E z increases 
from zero, the two peaks move oppositely in the TI region 
and then both move to higher frequency in the BI region. 
The plots shown here are for the total longitudinal con- 
ductivity. For the longitudinal conductivity about one K 
point, the same frequency dependence remains and one 
merely has half the magnitude, although the lower fre- 
quency segment would represent the response of charges 
with one particular spin orientation and the larger en- 
ergies are associated with the superposition of the re- 
sponse of both spin orientations. It is clear that from the 
frequency dependence, one could identify the spin orbit 
energy gap value by tracking the behavior with perpen- 
dicular electric field. The extrapolation of the minimum 
gap to zero as a function of E z identifies the critical elec- 
tric field that produces the VSPM, and the identification 
of the two regimes, TI versus BI, is found by looking for 
co-moving or counter-moving peaks. 

Of course, some things change in the doped case. Pri- 
marily, the lowest energy for absorption at a given elec- 
tric field strength is set by 2/i or A, whichever is greater. 
This is most easily illustrated with Figs. [5] and [6J where 
twice the chemical potential is taken to be finite but less 
than Ago in the former case and greater than Ago in the 
latter one. The 2/i Pauli-blocking cutoff is well-known 
in graphene where vertical transitions from an occupied 
state to an unoccupied state cannot occur in the band 
structure for photon absorption frequency < 2\/i\. In 
Figs. [S] and HO we see that the same pattern from the 
fx = case is occurring with the exception that there is a 
lower energy cutoff based on the inset figures illustrating 
the progression of the energies of the absorption edges 
with varying A z (following the arrows in the inset). For 
/i/Aso = 0.25, 2/i = 0.5Ago and the lower absorption 
edge decreases until f2 = 2/i where it remains stationary 
until it starts to increase in frequency at |A Z — Ago > 2/i. 
Note that an intraband Drude conductivity also appears 
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FIG. 6. (Color online) The real part of the dynamical con- 
ductivity Re<j xx (Q) /ao for /x/Aso = 0.75 and varying (a) 
A z /A S o < 1 and (b) A z /A S o > 1. Inset: A plot of the 
minimum and maximum band gaps, reproduced from Fig. 
showing the evolution of the absorption peaks with A z . For 
small Az there is a lower bound of 2/i before one peak emerges 
and then eventually a second peak. 




FIG. 7. (Color online) The Drude weight as a function of 
A z /A S o for varying fi/Aso- 



becomes strong enough to increase the minimum gap 
value further, the Fermi level will once again fall in the 
gap rather than in the spin split band and the Drude 
intraband absorption will disappear (as it does here at 
A z = 1.5 Ago). Thus, for 2/i < Ago, the Drude weight 
appears within the region where the minimum gap drops 
below 2fj, (as shown in the inset of Fig. [5]). For 2/i > Ago, 
the Drude weight is always present until the minimum 
gap overcomes 2/i (refer to the inset of Fig. EJ. The kink 
in the curve occurs when 2/i is equal to the maximum 
gap. Moreover, if the Fermi level is between the spin 
split bands, the Drude weight will show an increase in 
the TI region as A z is increased and then decrease again 
as the system is tuned into the BI region. 



IV. CIRCULAR POLARIZATION 



during the range of A z where the lowest absorption edge 
remains at 2/x, i.e. 2/i > |A Z — Ago I or 2/i > A m ; n . At 
this point the chemical potential is in the band and conse- 
quently intraband absorption processes can occur. This 
leads to an interesting observation that varying electric 
field at finite doping could provide a mechanism to switch 
on and off the Drude response at low frequency due to 
the Fermi level falling outside and inside the band gap as 
tuned by A z . This is illustrated in Fig. [7] where we plot 
the Drude weight (the amount of spectral weight under 
the delta function of Rea xx (tt) / ao for f2 = + ). Exam- 
ining the black curve with /x/Ago = 0.25, 2/i is less than 
the minimum gap A m j n = |Ago — A z for A 2 /Ago < 0.5 
at which point only interband transitions can occur. For 
increased A z , however, the minimum gap decreases and 
the Fermi level is in the band and an intraband Drude 
absorption occurs. The maximum occurs at the point at 
which the VSPM is reached and then the Drude weight 
decreases as the minimum gap once again increases with 
A z in the BI regime. When the electric field strength 



Circularly polarized light has been suggested as 
a mechanism for seeing the valley-spin coupling in 
silicone^ and in group Vl-dichalcogenide o 17 i 18 . In light 
of this suggestion, we discuss the full frequency depen- 
dence of the optical response of silicene to circularly po- 
larized light. Our discussion is somewhat different from 
that of other references^— where the focus has been 
on the matrix elements evaluated in momentum space 
for specific frequencies. Here we look at the broadband 
response. For finite frequency optical response in the 
circular polarization basis, the conductivity is given as 
a xx (n)±ia xy (n) for righthanded (+) and lefthanded (— ) 
circular polarization. 31 The real part of this function pro- 
vides the absorptive part of the optical response: 

a±(n) = Rea xx (n) T hntr xy (Sl), (17) 

which can be obtained from our previous formulas. This 
form refers to the total absorptive response to circularly 
polarized light which is written in terms of the absorptive 
parts of the longitudinal and transverse conductivities 
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FIG. 8. (Color online) The optical response to circularly 
polarized light for a TI with spin-split bands. The re- 
sponse is shown separately for each K and K' points for 
righthanded polarization, er+(fi)/<7o, an d lefthanded polar- 
ization, CT_(fi)/<7o- Here, A z = 0.75Aso and fi — 0. The 
blue shading is for the up-spin band, the red for the down 
spins, and the total response in one valley is the envelope 
of these curves following the black curve (sum of the red and 
blue curves) at higher frequency. The insets show the primary 
transitions occurring at the peak absorption frequency. 

which have been summed over valley and spin indices. 
It will facilitate our discussion to show the components 
which are summed over spin but plotted separately for 
each valley index, ie., 

4 (ft) = £[Re<f (0) T Im<(ft)], (18) 

and we decompose further to show the individual com- 
ponents based on valley and spin (the quantity in the 
square brackets of Eq. (fT8|l ). 

To begin with, we show in Figs. [5] and |H] plots to illus- 
trate the essence of our results. We have chosen for these 
plots to use the analytical formulas in order to simplify 
and sharpen the structure and to allow us to make com- 
parison with other figures later on for a discussion about 
the effect of broadening which we have not had up till 
now. In Fig. [51 we show the case of A z /Aso = 0.75 which 
is the TI regime and in Fig. [3] we take A z /Aso = 1-25 
which is in the BI regime. In the four frames, we plot 
the response to circularly polarized light for each com- 
bination of K and K' point and +/— polarization. The 
blue area indicates the spin up component and the red 
region, the spin down (as indicated by the arrows). 

Considering first Fig. [SJ it is quite apparent that with 
positive polarization and a frequency associated with the 
minimum gap fti = Ago — A z | = 0.25Ago, the optical 
response is dominated by charge carriers with down spins 
residing in the K' valley, whereas, for negative polariza- 
tion, at the same energy, the response is due to charge 



FIG. 9. (Color online) The optical response to circularly po- 
larized light for a BI with spin-split bands, at the K and K' 
points for righthanded polarization, a+(tt)/ao and lefthanded 
polarization, er_(n)/<ro. Here, A z = 1.25Aso and n — 0. The 
notation is the same as for Fig. [8] 



carriers with up spins in the K valley. Hence, circu- 
larly polarized light provides a means to populate valley- 
specific carriers with a particular spin orientation. The 
frequency required is that corresponding to the minimum 
gap in the band structure and indicated by an arrow on 
the band structure insets of the upper right frame and 
lower left frame of Fig. |5J Such an idea has been dis- 
cussed previously for silicene by Ezawai^ and by others 
for the dichalcogenides i 1 7 ' 18 It was also argued in those 
works that using a photon energy associated with the 
maximum gap would connect the two outer bands, thus 
allowing for absorption at that frequency which would be 
associated with the other spin orientation for each valley. 
That is, for 0,2 = Ago + A z = 1.75 Ago shown here, the 
dominate response is to see the down spin charge carriers 
in the K valley and the up spin ones for the K 1 point. 
However, as is clear from these four frames, the circularly 
polarized response at SI2 will not be purely of one spin 
orientation, nor when the response of the two valleys are 
added for a particular polarization, will it be the response 
of just one valley at that frequency. As shown in the in- 
set band structure pictures for the upper left and lower 
right plots, even though 0,% will connect the two outer 
bands at k = to give a peak in the absorption, it will 
also connect states from the inner two bands at finite k, 
and the total absorption is the sum of the two (the black 
curve). Thus, the lowest frequency £!i is robust for pro- 
ducing valley-spin polarized charge carriers but the upper 
frequencies will be contaminated and not result in a pure 
spin- valley polarized response. So we conclude that using 
photons with energy Oi will produce (if ,t) and (K',\.) 
in this case, given negative and positive circular polar- 
ization, respectively. The question might then be asked 
as to whether one could produce the other combinations 
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in an uncontaminated way. If one returns to the bottom 
frame of Fig. [3J it can be noted that if A z oc E z changes 
sign, this reverses the spin labels on the bands at the K 
valley (indicated by the dashed and solid lines) and like- 
wise at the K' valley and so, at least from a theoretical 
point of view, (K,\.) and {K' ,f) could also be excited at 
fii with oppositely polarized light if one could reverse the 
electric field without changing anything else. 

In the BI regime, shown in Fig. [9] the essential dif- 
ference is that the lower band spin response is switched 
between the two polarizations for a fixed K point and so 
a double peak structure appears in the positive polariza- 
tion response of the K valley and very little occurs in the 
negative polarization at the same point (with the results 
reversed at the K' valley) . Such behavior was also noted 
by Ezawa and is due to band inversion in the BI relative 
to the Tliiii Thus, one finds that once again, the lowest 
energy defined as Oi as before, will reflect a single spin 
state, but now associated with the opposite polarization 
from the case of the TI. For the total response to polar- 
ized light as a function of frequency, in the case of a BI, 
cr+ probes mainly the K point and er_, the K' one, and 
the two peaks in the frequency dependence are associated 
with different spin orientation. Whereas for the TI, <r+ 
excites primarily down spins and cr_, up spins, and it is 
two peaks in the frequency dependence that resolve the 
valley-dependence. 

Now we comment on the variation of these curves with 
A z . Given that the K 1 valley is always the same but re- 
versed with respect to the K valley, as we have discussed, 
we will only show the K valley results in the following. 
In Figs. [TO] and [TTJ we show the response for [i = for 
circularly polarized light a + (fl), <7_(i7), and the imagi- 
nary part of the Hall conductivity Ima xy (Q) as we vary 
A z . Note that -Ima xy (fl) = [<7 + (fi) - <r_(0)]/2 and 
is therefore a measure of circular dichroism. These fig- 
ures are only for the one K point. For the K' point, the 
curves are the same except reversed in polarization label 
and spin orientation information (see Figs. [5] and [SJ and 
reversed in sign in the case of — lma xy (il). The cases for 
the TI and VSPM are shown in Fig. [TO] and the BI, in 
Fig. [11] The first thing to note is that compared to Figs. [5] 
and [9] the peaks are suppressed and broadened due the 
impurity scattering included in these curves, Note, for 
instance, that the lowest peaks in Fig. [TOT b) are not at 
a value of 1 as seen in Fig. [5] but are reduced by a sig- 
nificant amount. This is due to the broadening and the 
reduction is stronger when the energy scale of the peak 
edge is smaller and therefore closer to the energy scale 
of the broadening parameter. As a result, these curves 
present a more realistic sequence in behavior with varia- 
tion with E z . While we do not identify the separate spin 
components as in Figs. [S] and O one can see that the over- 
all envelope of the curves in Figs. [5] and [3] (i.e., the sum of 
blue and red to give the black curve) is the same here and 
so one may use Figs. [5] and H] as a guide to identify the 
spin orientations associated with the peak structures in 
Figs. [TOI and [Tn Once again for the TI case, the pattern 




FIG. 10. (Color online) The optical response to circularly 
polarized light at the K point in the TI phase with /i = 
(a) Conductivity for righthanded polarization, a+(Q,)/ao. (b) 
Conductivity for lefthanded polarization, <j-(Q,)/ao. (c) The 
transverse Hall conductivity — lma xy (0,) which is also (o- — 
a+)/2. The peaks in (b) are associated with spin up charge 
carriers and those in (a) are predominantly due to spin down 
charge carriers. 

remains that the peak seen in the lefthanded polariza- 
tion cr_ is primarily due to up spins and that seen in the 
righthanded polarization er + is due to down spins super- 
imposed on a background of absorption due to transitions 
between the up-spin bands. In the BI case, there is only a 
response in one of the polarizations (righthanded in this 
case), however the lower peak is spin-selected for up spins 
while the upper absorption edge is again an admixture 
of the response of two spin species, dominated by down 
spins. The lefthanded polarization shows virtually no re- 
sponse by comparison. The interesting feature of these 
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FIG. 11. (Color online) As for Fig. [TT] but for the BI phase, 
(a) a + (Q)/ao. (b) cr-(fi)/W (c) -Imtr^ (fi) = (a- -cr+)/2. 
The peaks at lower energy are associated with spin up charge 
carriers and those at higher with predominantly spin down 



curves is that the absorption edges move with electric 
field tuning as seen before in the case of the longitudi- 
nal conductivity. As a result, the frequency at which the 
light can excite carriers with a certain spin orientation in 
a particular valley can be tuned (i.e., the f^i discussed 
before depends on A 2 oc E z ). Following the discussion 
in our previous section for the longitudinal conductiv- 
ity, there will be movement of these peaks in frequency 
according to the minimum and maximum gap and they 
will be suppressed by a 2fi cutoff when 2[i is greater than 
the gap energy. Such behavior will allow for additional 
flexibility in tuning the frequency at which valley-spin 
polarized carriers are excited, and the curves shown here 



illustrate that the response may not always be 100%. 



If one takes the difference between these two polariza- 
tions, i.e., (7+ — <r5 to obtain — 2Ini(7,^(Sl), a change in 
sign is seen for the TI case at the energy of fi 2 = A z +Aso 
but not for the case of the BI, where it is merely addi- 
tive. As we have seen, Q2 is the point where the dom- 
inant response transitions from one spin species to the 
other. Of course, due to the sign change from one val- 
ley to another in Imcr| y , the total response for Im(7 xy (f2) 
is identically zero. We note at this point that the to- 
tal <J±(Q), summed over valleys and spins, will be the 
same function as Rea xx (£l), but that the character of the 
absorptive response for circular polarization of light will 
be associated uniquely with a particular spin and valley 
index, particularly at the lowest absorption edge. Conse- 
quently, keeping this in mind, we do not show results for 
finite doping but use the longitudinal case for reference. 
The primary effect of finite doping where 2fi > fii is to 
eliminate absorption at lower frequencies and induce a 
Drude response. However, the Drude response will not 
show spin-valley polarization under circularly polarized 
light. 



Finally we discuss the VSPM. This is shown as the 
solid black curve in Fig. [TU] for A z /Aso = 1- The low 
frequency response of the Dirac-type dispersion is the flat 
background seen in (a) and (b) associated with up spins 
of the K valley. The magnitude is oo/4 in both cases. 
The result at K' is reversed such the flat background 
is due to down spins. As the two K points must be 
added for the total a±, there is an equal admixture of 
spin up and spin down, regardless of the sign of circular 
polarization and so one would need a mechanis m 32 ' 33 that 
would allow for sampling one valley only in order to see 
a single spin state associated with the linear bands. The 
only variation in response occurs at f2 > 2Ago when the 
upper band due to the opposite spin species enters and 
then a peak appears in a+ for down spins, superimposed 
on the flat background due to up spins. In — lma^ y (fl) at 
the one valley, there is only one peak as opposed to two 
for the BI or two with opposite sign for the TI with finite 
E z . As a final comment, we wish to note that for E z = 0, 
the bands are not spin-split and as a result, there is no 
circular dichroism or valley-spin polarization. 



It is important to reiterate that while we show valley- 
separated curves here, the total response of a particular 
polarization is the sum over the two valleys and hence, 
the result will simply look like the curves for Rea xx (£l) 
seen in previous figures. The point here is that the low- 
est peak in those curves is actually associated with a 
particular spin and valley when looked at with circular 
polarization. With other mechanisms 32,33 for separating 
out the valleys, the extra information shown here may be 
accessible. 
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V. SUMMARY 

Wc have calculated the optical response of silicene and 
related materials, such as germanene, which are predicted 
to exhibit a tunable band gap due to an applied perpen- 
dicular electric field. Due to the interplay of spin orbit 
coupling and the electric field strength, the bands dis- 
play spin- valley coupling. Tuning of E z , allows for rich 
behavior varying from a topological insulator to a band 
insulator with a valley spin-polarized metal at a criti- 
cal value in between. We find that in the longitudinal 
dynamic conductivity, two peaks marking interband ab- 
sorption edges will shift oppositely in frequency with in- 
creasing E z in the TI phase and in the BI phase, will 
move in tandem to higher frequency, separated by 2 Ago- 
The VSPM exhibits a flat graphene-like conductivity at 
lower frequencies changing to a single peak at 2Aso due 
to a second spin-split band with a band gap of 2 Ago- 
We suggest that seeing some of these features or track- 
ing the variation of the peaks with E z should allow for 
the identification of these phases and the determination 
of parameters such as Aso and the critical electric field 
that produces the VSPM. We have also noted the tun- 
able band gap of this system could allow for a switch- 
ing on and off of the Drude or DC response which may 
be of interest for technological applications. With Aso 
predicted to be about 1.5meV in silicene and 25meV in 
germanene, this spectroscopy should fall in the realm of 
the THz to far-infrared. These spectroscopic ranges have 
successfully probed the dynamic conductivity of a va- 
riety of materials, for example, the first observation of 
the far-infrared conductivity of graphene^ and the 1.4 
meV energy gap in superconducting Pb measured most 
recently in the THz range. 35 Consequently, this should be 
quite feasible once appropriate samples are developed. 

Another interesting aspect of this study is the conse- 
quences for circularly polarized light. We have shown 
that using circularly polarized light will resolve the spin 



and valley degrees of freedom. For instance in the TI with 
a perpendicular electric field yielding spin-split bands, 
the spin up charge carriers at the K valley can be re- 
solved by lefthanded circularly polarized light and the 
spin down ones at K' with righthanded polarization. 
Moreover, righthandcd/lefthanded polarization probes 
primarily spin up/down carriers as a function of fre- 
quency. Tuning the electric field into the BI phase causes 
a band inversion such that the A"- valley spin-up carriers 
are now resolved with righthanded polarization and spin- 
down ones at K 1 with the lefthanded case, and circularly 
polarized light resolves primarily one valley as a function 
of frequency. A reversal of the perpendicular field di- 
rection could also reverse the spin- valley polarization for 
these two cases. The range of frequency for obtaining a 
single spin orientation in the response is largely restricted 
to the region about lowest frequency peak in the inter- 
band absorption. Doping above the minimum gap edge 
will simply shift the spectral weight to an intraband or 
Drude component which will not be spin-valley resolved. 
At present, no experiments of this type exist on silicene, 
but a demonstration of tuning the band gap from TI to 
BI in such systems would be of great interest. Other 
materials, namely the group VI dichalcogenides, are in 
the BI limit and predictions for the dynamic conductiv- 
ity are being examined for M0S2 by others.— Given the 
potential for further developments with 2D crystals, we 
anticipate that the optical properties and band gap tun- 
ing will play an important role in the investigation of new 
materials and the development of new technologies. 
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